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We present a spin foam model in which the fundamental "bubble amplitudes" (the analog of 
the one-loop corrections in quantum field theory) are finite as the cutoff is removed. The model is 
a natural variant of the field theoretical formulation of the Barrett-Crane model. As the last, the 
model is a quantum BF theory plus an implementation of the constraint that reduces BF theory 
to general relativity. We prove that the fundamental bubble amplitudes are finite by constructing 
an upper bound, using certain inequalities satisfied by the Wigner (3n)j'-symbols, which we derive 
in the paper. Finally, we present arguments in support of the conjecture that the bubble diagrams 
of the model are finite at all orders. 
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O ' I. INTRODUCTION 

<n : 

To unlock the puzzle of quantum gravity, and understand quantum spacetime, we have to learn how to treat 
general relativistic quantum field theories in a background independent fashion jlj . A surprising number of research 
directions aimed at exploring background independent quantum field theory have recently been converging towards 
\^ , the spin foam formalism 

' A spin foam model can be seen as a rigorous implementation of the Wheeler-Misner-Hawking sum over 
geometries formulation of quantum gravity. The foam-like geometries summed over are spin foams, or colored 2- 
complexes. A 2-complex J is a (combinatorial) set of elements called "vertices" v, "edges" e and "faces" /, and a 
boundary relation among these, such that an edge is bounded by two vertices, and a face is bounded by a cyclic 
sequence of contiguous edges (edges sharing a vertex). A spin foam is a 2-complex plus a "coloring" N, that is an 
assignment of an irreducible representation Nf of a given group G to each face / and of an intertwiner i e to each 
\^ , edge e. The model is defined by the partition function: 

o: z = j2 m(j)Y1 U A *f Umn), (i) 

O ■ j n feJ v£j 

O i 

where Ajv is the dimension of the representation N, A V (N) is an amplitude associated to vertices: a given function 
of the colorings of the faces and edges adjacent to the vertex; and N(J) is a normalization factor for each 2-complex. 
(3J[v Let us mention a few of the research directions that, following very different paths, have converged to models of 
this kind. One of the oldest models of this type is the Ponzano and Regge formulation of 2+1 gravity pL0| . The 2- 
. £^ ° complexes in this case are the 2-skeleta of the dual of Regge triangulations, and the vertex amplitude turns out to be 
simply a Wigner SU(2) 6j symbol. In loop quantum gravity ]TT[ |, spin foams emerge as histories (in coordinate time) 
of quantum states of the geometry that is, histories of spin networks ||Jl|. In this case, the vertex amplitude 
A v {c) is given by the matrix elements of the Hamiltonian constraint. In covariant lattice approaches, the sum 
over colors corresponds to the integration over group elements associated to links, expressed in a ("Fourier") mode 
expansion over the group. In this case the vertex amplitude A v (c) is a discretized version of (the exponential of) the 
Lagrangian density |lj,|l5|. In topological field theories, the vertex function is a natural object in the representation 
theory of the group G, satisfying a set of identities that assure triangulation independenc e [|l7|-p0 |. Finally, in the 
modifications of topological quantum field theories related to quantum general relativity [Mlllnfl , the topological 
field theory vertex amplitude is altered in order to incorporate a quantum version of the constraints that reduces the 
BF topological field theory |^l|,^| to general relativity. 

Spin foams are very much analogous to Feynman diagrams. In references [ ^3| , ^4| , indeed, it is shown that the 
sum (Q) can be obtained as a Feynman expansion of a field theory over a group manifold, whose interaction terms 
determine the vertex amplitude (see also J25[). In this context, spacetime emerges from a Feynman expansion, as in 
the old 2d quantum gravity matrix models, or zero-dimensional string theory p3| . 

As in the Feynman expansions of a conventional field theory, one expects divergences to appear in the sum. There 
are two types of potential divergences: the ones associated to the sum over the colors of a fixed 2-complex, and the 
ones associated with the sum over 2-complexes. Here we consider only divergences of the first kind. In conventional 
Feynman diagrammatics, divergences are originated by integrating over the momenta circulating along closed loops, 
because momentum conservation at the vertices relates the momenta of adjacent propagators. In a spin foams model, 
the sum is over representations associated to faces, and a constraint analogous to momentum conservation is provided 
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by the requirement of the existence of non trivial intertwiners on the edges. Consequently, divergences are associated 
not to loops, as in quantum field theory, but rather to bubbles: collections of faces forming a closed surface ^7j. A 
way to control bubble divergences is to replace the group G with a g-deformed group, choosing q such that q n = 1, 
and to sum only over the finite set of proper representations of the quantum group. This is done, for instance, in 
the Turaev-Viro |^8| finite version of the Ponzano-Regge theory. The parameter q plays the role of a cutoff, and the 
physical theory is recovered by appropriately taking the q — ► 1 limit. 

The significance of the bubble divergences depends on whether the model is topological. A topological field 
theory is a diffeomorphism invariant theory which does not have local degrees of freedom, but only global ones. 
General relativity in 3 spacetime dimensions and BF theory are topological theories. On the other hand, a non- 
topological diffeomorphism-invariant theory is a theory, such as general relativity in 4d, which is generally covariant 
but, nevertheless, has local degrees of freedom (waves). In the context of spin foam models, the fact that a theory is 
topological is reflected in the fact that the amplitude of a fixed 2-complex in ([j]) is independent from the 2-complex 
(triangulation independence ones the manifold, i.e., the topology is fixed). More precisely, in general, in a topological 
model, bubble amplitudes diverges in the q — ► 1 limit, but topological invariance implies that the sum over colorings 
is the same as the one of a 2-complex in which the bubble has been removed, up to a divergent overall factor that 
depends only on the cutoff. In this case, therefore, divergent diagrams do not provide any additional information. 
The fact that a topological theory does not have local degrees of freedom is reflected in the triviality of all its 
"radiative corrections" (bubble diagrams). The consequence is that the sum over 2-complexes is trivial, and can be 
dropped, thus dropping all 2-complexes with bubbles. 

The situation is different in the non-topological context, and in particular, for quantum general relativity. In this 
case, the "radiative corrections" , that is, the bubble amplitudes, are the ones that carry the information about the 
quantum behavior of the local degrees of freedom of the theory. In this context the sum over 2-complexes is necessary 
in order to capture all the degrees of freedom of the theory, and bubble amplitudes are physically important. 

Examples of non-topological models are the covariant expansion in coordinate time obtained from loop quantum 
gravity, the Iwasaki jl(| and the Reisenberger [pi] models, and the Barrett-Crane model. The Barrett-Crane model 
is a non-topological modification of a topological model: 4d 5*0(4) BF theory, or the TOCY (Turaev-Ooguri- 
Barrett-Yetter) model |l7|-p0|]. It is well known that 4d 5*0(4) BF theory has a peculiar relation with 4d Euclidean 
general relativity: general relativity can be seen as a 4d 50(4) BF theory plus a constraint, which has an intriguing 
geometrical interpretation. The Barrett-Crane model is a modification of the TOCY model in which this constraint 
is implemented in the quantum theory |29|j3C| l. As a consequence of the constraints, topological invariance is lost and 
the model acquires local degrees of freedom. Radiative corrections associated to bubble diagrams carry non trivial 
physical information, but, unfortunately, diverge |27f . In order to be able to extract physical information from this 
model it is necessary to deal with these divergences. 

In a companion paper [p7| , we have begun a general study of the bubble amplitudes and the possible ways of 
renormalizing away their infinities (see also p4|). In particular, we have studied the version of the Barrett-Crane 
model developed in ]25| ] , in which the sum over complexes is explicitly implemented by obtaining the spin foam model 
from a field theory over a group manifold. In the course of this analysis, we have stumbled across a remarkable simple 
modification of the action of the Barrett-Crane model which leads to finite fundamental bubble amplitudes.^] We 
present this model in the present paper. We also argue that all "radiative corrections" in the model are likely to be 
finite. The model we present shares with the original Barrett-Crane model the feature that makes it an intriguing 
candidate for a quantum theory of (Euclidean) 4d general relativity: that is, it is another implementation of the 
constraint that reduces BF to GR. As a result, we obtain a theory, formally related to Euclidean quantum general 
relativity, presumably finite at all orders in the expansion over 2-complexes. 

The paper is organized as follows. In section II, we recall the field theory formulation of spin foam models by 
discussing the TOCY model ]l7|— pO| , and the Barrett-Crane model; then we introduce the new model. In section III, 
we study the model in configuration space. We define and compute the amplitudes corresponding to the fundamental 
bubble diagrams. We prove finiteness of the 1-bubble diagram. In section IV, we study the model in momentum 
space where we show that the interaction vertex of the theory is the Barrett-Crane vertex. We also give another proof 
of the finiteness of the 5-bubble amplitude. In the appendix we review some known facts about the representation 
theory of 50(4), and we prove certain inequalities satisfied by the Barrett-Crane vertex amplitude, the 6j, and the 
15j-symbols. We use these inequalities to show the finiteness of the fundamental bubble diagrams. They generalize 
to the 3Aj-symbols. 



1 The precise meaning of the fundamental bubble diagrams is explained below. They represent the simplest potentially 
divergent graphs in analogy to the one loop corrections in standard QFT. 
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II. DEFINITION OF THE MODEL 



We begin by recalling the formulation of quantum BF theory, and of the Barrett-Crane Euclidean quantum gravity, 
as field theories on a group manifold. These were first developed by Ooguri in |l8) and by DePietri et al in 
Then we introduce the new model. 

A. TOCY model as a QFT over a group manifold 



A spin foam model can be cast in the form of a field theory over a group manifold [£3 24 1. The simplest of these 
theories is given by the TOCY Jl7]-|o| topological model, corresponding to the quantization of 50(4) 4-dimensional 
BF theory. We begin by describing this model, which allows us to introduce definitions relevant in the sequel. The 
model is defined by the action. 

S[4>] = ^ J dgi . . .dg A 4> 2 {91,92,93, 9a) +^ j dgi . . . dg w 0(31,32,53,34) (2) 
0(34,35,56,37) 0(37,33,58,39) 0(59,56,32,3io) 0(3io,38,35,3i)- 

Here gi e 5*0(4) and the field is a function over 50(4) 4 . All the integrals in this paper are in the normalized Haar 
measure. The field 4> is required to be invariant under any permutation of its arguments; that is, </>(gi, 32, 53, 54) = 
0(3o-(i),3ct(2),3<t(3),3<t(4)), where a is any permutation of four elements; and under the simultaneous right action of 
50(4) on its four arguments: 

0(51,52,53,54) = 0(515,525,535,345) V.g G 50(4). (3) 
Let us introduce some simplification in the notation. We write 0(51,52,53,54) as 0(5i), and we write the (j^) as 



S'H = / dg l tt(g l )} 2 + ± I d5,[0(5,)] 5 , (4) 



where, notice, the fifth power has to be interpreted as in_ 

Instead of requiring that the field satisfies property 
</> and project it on the space of the fields satisfying (|3 
theory by the action 



j), we can also define the theory in terms of a generic field 
by integrating over the group. That is, we can define the 



5[0] = J d 9l [P s 0(5*)] 2 + ^ j d 9i [P s 0(5,)] 5 , (5) 
where the field is now arbitrary (except for the permutation symmetry), and the operator P g is defined by 

P g 4>{gi) = I ' d 1 0(5*7), (6) 



where 7 £ 50(4). 

Let us write the Feynman rules of the theory in coordinate space. To this aim, we write the action as 

£[0] = \ J dg l d~g l + ^ J dg^ V(gij) 4>{gij) 0(3 2j ) 0(53j) 0(34j) 0(35i) (7) 

in the last integral, i ^ j. Clearly, 4>{9ij) = 0(3i2, 3i3, 3i4, 315), and so on. 
The kinetic operator IC(gi,gi) can be written as 

>C(9i,9i) - E I d "l IPfo^-J)). (8) 

a J i=l 

The propagator is the inverse of the kinetic operator, in the space of the gauge invariant fields. The operator K, is a 
projector (i.e., JC 2 = K): its inverse in the subspace of gauge invariant fields corresponds to itself. The propagator 
of the theory is then simply 
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7 > (j9i,9i) = fc(9i,gii)- (9) 

The propagator is formed by 4 delta functions (plus the symmetrization and the integration over the group). Its 
structure can therefore be represented graphically as in Fig. flL in which every line represents one of the delta functions 
in §). 



FIG. 1. The structure of the propagator. 
The potential term can be written as 

V(9ij) = ^ [dfo J] S( gij Z?- 1 /?^ 1 ), (10) 
i<j 

where Pi G 5*0(4). By introducing the notation pji = gj^gij, which satisfies the property pj^ = pji, and rearranging 
the arguments of the delta functions using the fact that 5{g\g2) — <5(32<?i), we can write the potential term as a 
function of ten variables only 

v(pij) =n/*n mpji^ 1 )- en) 

The structure of the vertex is represented in Fig. ||, where each line represents one of the delta functions that appear 
in the previous expression. 




From now on we will use the Greek letter (3 to denote the symmetrization integration variables associated to 
vertices, while we reserve the letter 7 for propagators. 

The Feynman diagrams of the theory are obtained by connecting the five-valent vertices with propagators. At the 
open ends of propagators and vertices there are the four group variables corresponding to the arguments of the field. 
For a fixed permutation a in each propagator, one can follow the sequence of delta functions with common arguments 
across vertices and propagators. On a closed graph, each such sequence must close. By associating a surface to each 
such sequence of propagators, we construct a 2-complex [ p3[ . Thus, by expanding in Feynman diagrams and in the 
sum over permutations in (ph, we obtain a sum over 2-complexes. In other words, a 2-complex is given by a certain 
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vertices-propagators topology plus a fixed choice of a permutation on each propagator. From now on, a "Feynman 
diagram" , or simply a diagram, will denote one such 2-complexes. 

According to the standard Feynman rules deduced from the form of the action, the amplitude of any diagram 
will be given as an integral over the internal coordinates -in this case SO (4) group elements- of the product of the 
corresponding propagators and vertices. Thus, the amplitude of any Feynman diagram is given by an integral of 
a product of delta functions over a (compact) group manifold. From (||) and ([Io|) we observe that the integrand 
defining this amplitude correspond to a product of delta functions related to each other according to Fig. [l] and |[ 

Working in configuration space, one can easily prove the topological invariance of the theory using the Feynman 
rules derived from (||) J27j|l^]. If one expands the action (g) in terms of irreducible representations using the Peter- 
Weyl theorem, then the amplitude of a Feynman diagram J is 

^( j )=En A ^n i5 ^)> ( i2 ) 

N f v 

where 15j(N v ) denotes the 15j-symbol constructed out of the 10 colors of the surfaces meeting at the vertex plus 
the five colors of the intertwiners corresponding to the five edges defining the vertex. The sum goes over all possible 
compatible colorings of faces and intertwiners in the 2-complex. 



B. The Barrett-Crane model in coordinate space 

Now we describe the Barrett-Crane model as a field theory over 50(4) 4 . Consider the fundamental representation 
of 50(4), defined on 5ft 4 , and pick a fixed direction t in 5ft 4 . Let H be the 5*0(3) subgroup of 50(4) that leaves t 
invariant. We define the projector Ph 

Ph4>(9i) = / dhi <j)(gihi), (13) 



where hi £ H. From now on, the letter h will always indicate an element in H = 5*0(3), while we use g, (3, 7, and p 
for elements in G — 50(4). The Barrett-Crane model is given by the action: 

5|>] = J d 9i [PgPh^gi)] 2 + ^ fdgi [P g Ph4>{9i)f , (14) 

where the notation is as in the previous section and the fields are assumed to be symmetric under permutations of 
their four arguments. In |23|| , it is shown that the Feynman expansion of (|l4|) gives the Barrett-Crane spin foam 
model, summed over 2-complexes. 

Notice that the action (|lj) is simply obtained by adding the Ph projection to the the action (||) of the topological 
TOCY. The projector Ph projects the field over the linear subspace of the fields that are constants on the orbits of 
H in G, that is, that satisfy 

<X5i: 52, 53, 54) = Hgihi, 32^2,53^3, g±K) VhiEH. (15) 

When expanding the field in modes, that is, on the irreducible representations of 50(4), this projection is equivalent to 
restricting the expansion to the representations in which there is a vector invariant under the subgroup H (because 
the projection projects on such invariant vectors). The representations in which such invariant vectors exist are 
precisely the "simple" representations, namely the ones in which the spin of the self dual sector is equal to the spin 
of the antiselfdual sector. In turn, the simple representations are the ones whose generators have equal selfdual and 
antiself dual components, and this equality, under identification of the 50(4) generator with the B field of BF theory 
is precisely the constraint that reduces BF theory to GR. Alternatively, this constraint allows one to identify the 
generators as bivectors defining elementary surfaces in 4d, and thus to interpret the coloring of a two-simplex as the 
choice of a (discretizcd) 4d geometry H]7],|3(J . 

Using equations (^), and ( |l3| ) it is straightforward to compute the kinetic operator and the interaction vertex of 
the Barrett-Crane model in coordinate space. The kinetic operator of the theory is 

^(ff»,5i) = J d^dh l ^5{g i hi'yh l gl 1 ). (16) 

i 

The vertex is 
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V(gij) = i / dfadhj Yl S(g ji 1 h lJ f3 i 1 /3 j h ji g ij ). 

i<j 



(17) 



Notice that the kinetic operator K. is not a projector anymore (i.e., K? ^ K,). As a consequence, the propagator 
V = Kr l docs not have a simple form in coordinate space. 

If we expands the action ( |l4| ) in terms of irreducible representations using the Peter- Weyl theorem, then the 
amplitude of a given Feynman diagram J is given by 



^ j )=En dim ( c /)n^)< 

C f V 



(18) 



where the sum is now only over simple representations of SO(A), and B(c v ) denotes Barrett-Crane vertex amplitude. 
This is given by the 1 5.7-symbol con structed out of the 10 colors of the surfaces meeting at the vertex using the 
Barrett-Crane intertwiners. See (A7) and ( |23||). 



C. The new model 



The idea at the basis of the new model is the same as in Barrett-Crane: to modify the BF action (||) by inserting 
the projector ( p"3| ) which implements the restriction of BF to GR. However, this time we insert the projection in the 
interaction term only, keeping the same kinetic term as in the BF theory. That is, we define the new model by the 
action 



S[4>] = J d 9i [P g 4>( 9i )] 2 + d 9l [P g P h <t>{9i)f 



(19) 



d 9i [ P g<i>{9i)Y 
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d 9i [p g p h PM 9i )Y 



(20) 



where Ph and P g are defined in (g), and (^3|) respectively. As for the Barrett-Crane model, if we drop Ph from 
the previous action we obtain the TOCY topological model of section (II). As we show below, the two forms of the 
action define the same theory, since the extra P g in the second expression can be always absorbed into the P g of 
some propagator when computing an amplitude. The second form of the action simplifies the analysis of the theory 
in momentum space. 

From (n3) the kinetic operator becomes 



K>{9i,9i)= J d^fdhi U%i75i 



(21) 



which corresponds to the projector into the space of gauge invariant fields [K? = 1C). In this space its inverse is itself 
and the propagator of the theory is simply V = JC as in the TOCY model. The vertex of the theory is 



vow. 



i 

5! 



dPtdfiidh^ Y[ S (9ji 1 Pi h ijPi 1 Pj h aPj l 9ij)- 



(22) 



The (3 and (3 integration variables correspond to the two projectors P g in the interaction. In the following we will 
show that the (3 integration is redundant when computing any amplitude, and therefore the two expressions in ( |l9| ) 
and (|2tj| ) define the same theory. The key feature of this model is that, as we will argue in the following, it does not 
contain divergences associated to the sum over colors for fundamental bubble diagrams. 



III. BUBBLE AMPLITUDES 



A. Fundamental bubbles and Pachner moves 



As mentioned in the introduction, and discussed in detail in Ref. |27 
in Feynman diagrams containing bubbles, that is, closed surfaces. Bub 



, in a spin foam model divergences arise 
jles are the spin foam analog of the loop 
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diagrams of standard QFT. To begin the analysis, we consider here only 2-complexes that are (the two skeleton of 
the) dual of regular triangulations. Therefore, we will talk equivalently about triangulations or 2-complexes. 

We define a fundamental bubble diagram as a bubble diagram obtained from an elementary diagram without 
bubbles by means of the basic 4d Pachner moves on the corresponding triangulation. These are the simplest diagrams 
presenting bubbles, and, in this sense, fundamental bubble diagrams are analogous to one-loop diagrams in standard 
QFT, representing the basic potentially divergent amplitudes in the model. 

In four dimensions there are 3 possible Pachner moves -the 1-5, the 2-4, and the 3-3 Pachner moves respectively- 
plus their inverses. Only the 1-5 and the 2-4 Pachner moves generate bubbles. In terms of a 4d triangulation, the 

1- 5 move is defined by the split of a 4-simplex into five 4-simplices. More precisely, we put a point p in the interior 
of the 4-simplex with vertices pi,i = 1, . . .5, we add the five segments (p,Pi), the ten triangles (j),pi,pj), the ten 
tetrahedra (p,Pi,Pj,Pk), and the five 4-simplices (p,Pi,Pj,Pk,Pi) (where i ^ j ^ k ^ I), 

In the 2-4 move, we replace the two 4-simplexes {a,Pi,P2iPz-,Pi)i and (b,Pi,P2)P3)P<i)> sharing the tetrahedron 
ipiyP2,P3,P4) with the four 4-simplices (a, b,pt,pj,pk) where i ^ j k = 1, . . .4. 

In terms of the 2-complex (the 2-skcleton of the dual of the triangulation) which represent Feynman diagrams of 
our field theory, this set of moves generate the two fundamental bubble diagrams. The 1-5 move in the dual picture 
is illustrated in Fig. [| We denote the diagram on the right as the 5-bubble diagram. The vertices of the picture are 
dual to the 4-simplexes of the triangulation; the edges of the picture are dual to the tetrahedra of the triangulation; 
the surfaces (in fact, here, all the triangles) are dual to the triangles of the triangulation. 

The amplitude of a closed diagram is a number. The amplitude of an open diagram, that is, a diagram with a 
boundary, is a function of the variables on the boundary, as for conventional QFT Feynman diagrams. The boundary 
of a 2-complex is given by a graph, where the nodes are generated by the intersections of the edges with the boundary, 
and the links are generated by the intersections of the surfaces with the boundary. To start with, the amplitude of 
the open diagram is a function of 4 group arguments per each external line. However, consider a surface of an open 

2- complex and the link ab of the boundary graph that bounds it. Let a and b be the nodes on the boundary graph 
that bound ab. The surface determines a sequence of delta functions that starts with one of the group elements in a, 
say g a and ends with one of the group elements in b, say gb- By integrating internal variables all these delta functions 
can be contracted to a single delta function of the form 5(g a ■ ■ ■ g^ 1 )- This is of course a function of g b ~ 1 g a . We can 
thus define the group variable p a b = g b ~ 1 g a ,: naturally associated to the link ab, and conclude that the amplitude of 
an open 2-complex is a function A(p a b) of one group element per each link of its boundary graph. In "momentum 
space", the amplitude of the diagram is a function over the possible colorings, in the sense of the spin networks, of 
the boundary graph. That is, if s is a spin network given by a coloring of the boundary graph, 

A(s) = J dp a b lp s (pab) A(Pab), (23) 

where V's(Paf)) is the spin network function p2| . In Fig. |^, the thin lines in the picture represent the boundary graph 
of the diagram, that is the intersection of the 2-complex with a 3-sphere that bounds it. This intersection is a graph 
on the 3-sphere. Notice that the boundary remains the same after the move is implemented. We can think of the 
diagram on the right as a radiative correction to the diagram on the left. The analogous picture is shown in Fig. ^ 
for the 2-4 Pachner move. We denote the diagram on the right of Fig. as the 1-bubble diagram. 




FIG. 3. The 1-5 Pachner move, on the right the fundamental 5-bubble Feynman diagram. 
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FIG. 4. The 2-4 Pachner move, on the right the fundamental 1-bubble Feynman diagram. 



A straightforward calculation shows that the amplitudes of these diagrams diverge in the cases of the TOCY and 
Barrett-Crane model j^jj. In the following section, we compute the amplitudes of the 5 and 1-bubbles. 



B. The 5-bubble amplitude 

We start by computing the amplitude A 5 corresponding to the diagram on the right of Fig. |[ We use the second 



form of the action (20) and show that the result is independent of this choice. The amplitude is a function of the 
group elements on the external propagators. It is obtained by writing the five vertices and the ten propagators in the 
diagram, and integrating on each propagator-vertex contraction. Notice that we are computing the amplitude of a 
fixed 2-complex, and therefore the permutation in (^]) is already fixed, and we do not have to sum over permutations. 

Let us introduce some notation. We label the five vertices in Fig. [| by means of an index a that takes values from 1 
to 5. The ten internal edges can then be labeled by unordered couples (ab) of indices. For notational convenience we 
label as (0a) the external edge emanating from the vertex a. Consider the a vertex. It is a function of twenty group 
elements, four per leg, naturally paired by being in the same delta function, or, equivalently, by relating to the same 
2-complex surface. Consider the leg (ab) of the vertex a. We denote the group elements on this leg which is paired 
with the leg (ac) as gi a b)( ac )- Notice that the first couple of indices refer to the leg on which the group sits, while 
the second refers to the leg to which it is paired. The four group elements on external lines are denoted as g(o a )(ab), 
while the group element on, say, the (ab) edge and paired to the external edge is denoted <7( a &)(a0)- Since there is 
one h integration variable corresponding to each of the configuration variables <?( a 6)(b c )j we label them according to 
the same rule. We denote the 7 integration variables in ( pl| ) associated to propagators connecting the vertex a with 
the vertex b as jab respectively. Given a vertex there is one (3 and one p per leg corresponding to the two P g 's 
in the interaction (see (pOj)). Therefore, the /3 integration variables in ( p2[ ) are denoted as f3 a (ab), where the first 
index denotes the vertex to which (3 belongs, while the second couple of indices denote the corresponding leg. The 
amplitude is then given by 

M(9a(ab)) — / [dg] n v (9(cb)(b d) ) II v{ 3i 

cd)(ce),g(cd)(de)), (24) 

b cd 

where [dg] denotes the integration over internal group elements (namely, g( a b)(bc) with a, b 7^ 0), while the indices 
c, d, e in the vertices and propagators can take the value as well. As argued in the last section, the 3(o a )(ab)' s enter 
this expression only in the combination p ab = .9^) (ab )5(o6)(a6)- 

We have to insert the values (|2^) of the vertex amplitude and the value (^,^) of the propagator (without the sum 
over permutations) in this expression. Explicitly, and with the correct index structure we need, these are 

V(5(c6)(6<0) = T\ jW]W][dh] \\ 5(ff^ )(cb) /36(c6)/«(c6)(6d)/3 6 ^6)A(6d) ft (6rf)(c6)/3 6 " ( b d) 3( c 6)(6d)) (25) 



(cb)<(bd) 



and 



'P(9{cd)(ce),9(cd)(de)) = I [dry] Yl S (9(cd)(ce)Jcdg (cdKde) ), 

e— 0,... ,5 7 e^c,e^d 



(26) 



where [df3], [df3], [dh], and [dj] denote integration over all the /3 a (ab), $a(ab)i ^(ab)(6c)i an d lab variables respectively. 
We recall that all h's are in 5*0(3), while all /3's, 7's, g's, and p's are in SO(i). Inserting ( |25| ) and ( p6| ) in ( ]2~i| ) we 
obtain a complicated multiple integral of delta functions, which we now analyze by breaking it into pieces. 

The key observation is that the delta functions appear in sequences, corresponding to the boundaries of the surfaces 
of the faces of the 2-complex, or to the connected lines obtained by replacing vertices and propagators in Fig. || with 
the vertex and propagators given in Fig. |l| and |2|. It is easy to see that there are two kinds of such sequences, 
corresponding to the two kind of surfaces in the 2-complex. These are illustrated in Fig. 5. 




FIG. 5. Two possible terms contributing to the A$ amplitude 



The terms of the first kind follow the pattern shown in the diagram on the left of Fig. 5; we denote them as 
wedges. The ten terms of the second kind have the form of the diagram shown on the right, and we denote them as 
faces. There are ten terms of the first kind and ten terms of the second kind. Consider the term of the first kind. 
Let (ab) be the link in this term. From the definitions of the vertex and propagator respectively we have that the 
corresponding sequence of delta functions is 

^(9(0a)(a&)/3 a ( a)^(0a)(afc)/3 a( L)' 3 a(afc)^(a6)(0a)/3 a (l lfc ).9( Q 1 b) ( 0a ))'5(5( af) )(0Q)7(a6)3( a h)(06)) 

^(ff(a6)(06)/3^ah)^(06)(a6)/3^fc)/ 3 KOb) fc (0^)(^)^KOb) 3 (W)(afc))' ( 27 ) 

where the first and the third delta functions come from the vertex a and the vertex b respectively, and the delta in 
the middle comes from one of the four deltas in the propagator ab connecting the two vertices. By integrating over 
the two variables that concatenate the delta functions, namely gr a b)(oa) an d ff(oh)(06)j we obtain the quantity 

Eab = KPb(ab) h (ab)(m)Pb(lb)Pb(Ob) h (Ob)(ab)^ b) Rba^ (28) 

We have rearranged the argument of the delta function using that ^(31.92) = <5(.92<7i), and we have defined p a b = 

9(0a){ab)9(0b)(ab)- 

The second kind of terms, illustrated on the right of Fig. Ill B| , come from the ten internal faces in the 5-bubble 
diagram. Consider the face bounded by the three vertices abc. The corresponding term gives the sequence of deltas 

S (9(ac)(ab)$a(ca) h (ac)(ab)P~( ca) P a (ab) h (ab)(ac)P^^ 
5 ^9(ab)(bc)P b (ab) h (ab)(bc)^ ( \ b) (i b ^ c) h ( ^ 

5 {9(bc){ac)P c( bc) h (bc)(ac)l3~^ c) l3 (c)ca h {a ^^^ (29) 

Again, we can integrate over the intermediate group elements, g( a c)(ab), 9(ab)(ac), g(ab)(bc), 5(6c)(o6), 9{bc){ac), and 
g(ac){bc) obtaining 

F abc ee <5(/3 a(afc )/i (a6 ) (QC )/3^ &) 7 (a6) /3 b(Qb) /i (a6)(; , c) /3^ b) Pb(bc)h( bc )(ab)$ b{ l c) 



Pit 



7(6 c )/3 c (6 c )^(6c)( a c)^ c (L)/ 3 c(c a )ft-(ac)(6c)/3 c( c a )7(ac)/5a(ca)ft-(ac)(afc)/3 a( c a) ), (30) 
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Notice that the 7 a b's associated to the internal propagators appear surrounded by /3's in the same way in E a b and in 
F a b c ] therefore we can reabsorb all the /3's by redefining the integration variables j a b (that is ^ a b — * $r a ) a bYabP(b)ab) 
using the right-left invariance of the Haar measure. This shows that the P g on the right in the interaction term of 
the action ( pp| ) is redundant and that the two forms (|l?J) and (|o|) of the action are equivalent. 
Using the definitions (53) and (^9|), we can rewrite the amplitude (153) as 



A 5 (p ab ) = J [dry] [dp] [dp] [dh] J] E ab ]J F abc . (31) 

ab abc 

We can eliminate the E a b's by integrating over the 7 a fc's. This gives 

lab = /3a(a6)Vfc)(Oa)/3~ ( * &) /3 a(0a) fr(Oa)(^ (32) 

We substitute these values into the remaining deltas; we define (see the emphasized term in equation (|29|)) 

Pab = Pa(ab)hab)(ac)h( a b)(0a)/3 a ^ b )P a (0a)h0a)(afy ( 33 ) 

We can redefine the /i( a b)(b c ) variables, by absorbing the /i( a 6)(ob) m to them. We substitute these expression in the 
F's, and the amplitude becomes 

MPab) = f W] [d$] [dh] \{8{p^ P t c 'pta). (34) 

J abc 

In order to show that A§(p a b) is finite, it is sufficient to show that the "vacuum" value of the amplitude is finite 
p7| . In momentum space this is the value of the amplitude when the colors of the boundary lines are zero. In 
configuration space this value is obtained by integrating the amplitude (|34| ) over p a b (Integration over the external 
group variables projects the amplitude to the trivial representation). We want thus to show that 

J 5 = J [dp] A 5 (p ab ) < oo. (35) 

The amplitude when the external colors are not all zero is given by the integral of A < ^' c \p) times the appropriate 
spin-network state, which is a bounded function over the group. The result of this integration is finite if ( J3q ) holds. 

If we insert (33) and (|34|) into (^H|), we see that several of the variables can be absorbed into redefinitions of 
integration variables. In particular, we can redefine p ab -> /3 o ( 0o )' l (0a)(a6)/5o(0a)' 9 « b A(0b)' l (06)(a6)^ 6 " ( 06) in ©i and 
obtain the following simpler expression 

[dp] [dp] [dh] JJ^(^ a(ai ,)/l(a6)(6c)/3^i 6) Po6/3 6 ( a 6)/l(c6)(a6) ) 3 b " ( ^ ) 
abc 

Pb^htabKb^Pb^PbcPcibcfiibc^P^bc)! 3 ^ ( 36 ) 



We can write this as 



where 



T 5 = J [dp] [dp] [dh] \{5{ptbKcKd ), (37) 

abc 



Pab - P(a)ab h {ab)(bc)P {a ) a bPabP(b)abh(bc)( a b)P (b )ab- ( 38 ) 

Now, in order to simplify this integral, we observe that the integrand in (^6|) is gauge invariant under the gauge 
transformation p a b — » p a PabPb~ X with p a £ S'0(4), since the p^s can be absorbed by redefining the /3's. The integral 
is therefore equal to the gauge fixed integral at a particular gauge, times the volume of the gauge orbit (which is unity 
since we are working with the normalized Haar measure). We fix the gauge by requiring pi2 = P13 = pn = pi$ = 1. 
Taking this into account, we can integrate away the variables p^-, where an index with hat takes the values 2 to 5 
only. The delta functions -associated to faces in which 2 of the 3 p's have been gauge fixed- fix the value of these 
variables to 
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P'bc ~ Pb(b§ h (bc){cl)P b{b& ^ 

/ 3 i(ai)^(ei)(i6)/ 3 i(6i)/ 3 e(ei)' l (i6)(ci)/ 3 c(ci)/ 3 c(Sc) ft '(ei)(Sc)/ 3 e(ge)' ( 39 ) 

Inserting this in (^8|), the external (3's simplify and two K can be collapsed into one by redefinition. There remain 
only the four delta functions associated to the faces abc. Each h variable turns out to be sandwiched between (3 and 
/3 _1 . Thus, integrating in such an h variable is indeed integrating over the subgroup f3H of 5*0(4), formed by the 
elements that leave (3t invariant. Taking all this into account, we conclude 

I5 = / df3 b{cd) / df)(oi>)(c<0 J! s (pltPb b cPcd ): (4°) 

J JP(ab)(cd)H abc=2,...,5 

with 



Pat = b(bc)(cd)§(ab)(bc)§(ab)(ac)t)(ca)(ab)t)(ac){bd)§(bc)(bd), (41) 

where t)( a fc)(b c ) S 0( a b)(bc)H (each t) corresponds to a different S'0(3) subgroup). 

As a check, imagine now that we are dealing with BF theory. To obtain BF theory, we have simply to drop 
the integrations over /i's. The resulting divergence is immediately (<5(1)) 4 , which is the correct divergent factor for 
the 1-5 pachner move in BF theory jl8[,^7|]. Does the integrations over the subgroup H absorb all the divergences? 
If it wasn't for the integration of the /3's, the answer would be no, because we would simply obtain a divergence 
proportional to the fourth power of the value on the origin of the delta on G/H. However, the combined integration 
over the ft's and the /3's is sufficient to absorb all the divergences. 

In fact, the integral we are studying is over a compact domain. Therefore divergences can only come from 
particular points where the integrand diverges. The point where the most serious divergences can appear is the 
origin f3( a b)( c d) — h(ab)(cd) = !• Let us study the behavior of the integral around this point. To this aim, it is 
sufficiently to study an arbitrarily small neighborhood of this point. On an arbitrarily small neighborhood of the 
origin, we can describe the group by means of its algebra. We can thus replace the group integral with an integral 
over the algebra. To the order we are interested, we can take the Lebesgue measure over the algebra, since the 
Jacobian between this measure and the measure induced by the Haar measure on the group goes smoothly to 1 on 
the origin. Group elements can be written as exponents of algebra elements, and products of group elements can be 
expanded to first order around the origin, giving commutators in the algebra. The four SO (4) delta functions give 
24 1-dimensional deltas. It is lengthy (because of the many variables), but completely straightforward, to see that 
the 24 resulting 1-dimensional delta functions in the previous expression are not redundant and therefore the value 
of the amplitude is finite. We will present the details of the calculation elsewhere. 

A simple proof of the finiteness of the previous amplitude will be given in terms of the mode expansion in momentum 
space. 



C. The 1-bubble amplitude 

We now consider the amplitude of the diagram on the right of Fig. f|. This is another divergent amplitude in BF 
theory and in the Barrett-Crane model. We now show that the amplitude is finite in the new model. The pattern 
for analyzing this amplitude is the same as in the previous case, but simpler. We only sketch here the key steps of 
the calculation, leaving the details to the reader. Integration over the six upper external wedges (see Fig. |j) implies, 
in analogy to the previous case, that 

lab = Pa(ab) h (ab) (Va)P a ( ab )P a (Qa) h (0a) (ah) a (la)P^Pb(Ob) h (0b) (a&) A,(0&) Pb(ab) h (ab) (Ob) P b {ab) ' ( 42 ) 

substituting this into the remaining deltas we obtain (schematically) 

A i(p) = JW][dh] n^(/3a(afc)^(BA)(OA)/3^ fc) /3 a(0a) ^(Oa)(Q6)Pa6^(06)(a6)/3 fc( o fc) 
abAB 

Pb(ab) h (BA)(OB)P b (l b )Pb(OB)h(0 B )( BA )PBAh( OA )( BA )P a ^ A) ) 

Yi S (Pab ^ PbT } Pea C) ) II 6 (( 3 a(Oa) h (Oa)(OA)PaAh { OA)(Oa)Pa(OA))^ ( 4 3) 
abc aA 
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Now we want to show that the "vacuum bubble" is finite, namely: 

X 4 = / [dp] Ai(p) < co. (44) 

Integrating over pab and p a A this reduces to 

J 4 = J [dp}[d0\[dh\Y[^pTpt: c) P^ c) )- (45) 

abc 

As in the previous case we can fix the gauge by means of the conditions pi2 — pn = pu = 1. We then integrate 
over the remaining p's which eliminate three of the four delta functions, and fix the value of the remaining p cc i (for 
c, d 5^ 1) to 

Pbc = A,( bc ) h (bc)(lb)Pb{bc) Pb(ab) h (ab) {ac)Pb{ab) @a(ab) h(ab)(bc) P~( ab ) 

Pa{ca)h{ca)(ab Pc(ca) h (ca)(cb)P c (l a) l3 c{hc) h( bc){ca) f3 c{ l c} . ( 46 ) 

Finally, reduces to an integral of single delta function. This integral is the integration of the distribution 



^(Si? 92, 93) defined in (A8). By lemma(l) in the appendix, the integral is finite. 



IV. SPIN FOAM FORMULATION 

Using Peter- Weyl theorem, we can analyze the new model in momentum space. This analysis explicitly connects 
our results to the spin foam formalism. In particular, we show that the vertex amplitude of the theory is essentially 
the same as the one in the Barrett-Crane model, up to a rescaling, that can be viewed as an edge amplitude. We also 
present here another proof that the amplitudes calculated in the previous section are finite, and finally we present 
arguments supporting the conjecture that amplitudes are finite at any order. 

According to the Peter- Weyl theorem, given an £ 2 [G] (Haar-square integrable) function 4>{g) over the group 
G = 50(4), we can expand it in terms of the matrices D^)(g) of the irreducible representations A of SO(A); that is, 

#S)=£ ^pD^ig). (47) 

A 

We begin by analyzing the kinetic term in ([H]) . Using @ and equation ( |A4| ) we obtain 

p^( 3l ,..., 54 )= £ ^■■■»4l:% ^ l)71 (9i)---^ 4)74 (.9 4 )E C7 ( A )7 1 ...7^(A)ft-A, (48) 

JV1...JV4 A 

See the appendix for notation and definitions. In order to have a more compact notation, we have dropped the indices 
N\ . . . N4 from the expression for the intertwiners and we have kept only the color A that labels them (namely we 
use Cy 1 1 ''j^ i A = C(A) 7l ... 74 ). We use Einstein summation convention over repeated indices. There is no difference 
between upper or lower indices and we write them in a way or the other for notational convenience. We can simplify 
the previous expansion by defining the new field components 

ai...o 4( oM \/3i.../3 4 



V ATi...AT 4 ,A — / a a \ — ' 



where Ajv denotes the dimension of the irreducible representation of order N, and the factor (A^ . . . Ajv 4 )^ has 
been chosen to simplify the expression of the interaction vertex computed below. This choice of field yields the mode 
expansion: 

P g 4>{g x , . . . ,g A ) = J2 • ■ ' A ^) § '"v :'V,a ( 9l ) . . . D^( 9i ) C(A) 7l ... 74 . (50) 

Ni...Ni,A 



Using (Al) and the orthonormality of the intertwiners, the kinetic term in (|19[) becomes 



K= E ^Xa^n^J^ ■ . . A N4 f 5^ . . . <W 4 . (51) 

Ni...N 4 ,A 
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We can directly read the propagator of the theory from this expression 



{/\ Nl . ../\ Ni r 



In order to write the potential term w e ne ed to expr ess P g PhP g 4> in terms of irreducible representations (see (19)). 
Starting with (Ea) and using equations flA4|) and (A5) we obtain 



P g P h P g c» = ]T (A Nl ...A N J $%\;XA D £ lhl (9i)---D^( 9 ±) (53) 



Ni...Ni,A 



C{N) 11 ... li C{N)^-^w Pl . . . w Pl w^ . . . w^C(A)^. 



/ j - \- /71---74 
N 



Applying equation (A6) and ( |A7| ) we obtain 

P g P h P g 4>= VA^.-.A^ ^V.VZ,A^ lhl (.9i)---^ 4)74 (.94) S 71 ... 74) (54) 

ATi...JV 4 ,A 

where the sum is now over simple representations only, and -B 7l ... 74 denotes the Barrett-Crane intertwiner. Using 
the previous equation the potential term in ( |l9| ) becomes 

' \ " \ " tf,a 1 a 2 a 3 a4 1 .1,04050607 .^07030309 ,1,090602010 ^.010080501 r> 

gj 2-^i 2-~i NiN 2 N 3 N i ,A 1 ^ NiN^NaNjM N 7 N 3 NfsN 9 ,K^ N 9 N & N 2 N 10l A4 ^ JVi Af 8 Ar 5 Afi ,A 5 D Ni_,... ,N 10 , \<JO) 
Ni...Nio A1...A5 

where Bn ± ,... ,n 10 corresponds to a 15j-symbol constructed with Barrett-Crane intertwiners which corresponds to the 
Barrett-Crane vertex amplitude M . Explicitly, 

n f,N 1 N 2 N a N 4 . tdN&NsNqN? r,N 7 N a N 8 N 9 r,N 9 N 6 N 2 N 10 riNmNaN^Nt /r fi \ 

^/Vi,... ,iVio ■ ^01020304 Q4Q5CK6CK7 07030809 090902010 Q10O8O5Q1 ■ v / 

Thus, the potential part of the action in the new model is given by (|55|) as in the Barrett-Crane model. Notice, 
however, that there is an extra sum over A in (^H|), absent in Barrett-Crane. The propagator (^2|) of the theory in 

momentum space is rescaled with respect to the Barrett-Crane propagator (P Ql)11 ... a4/J4 = 

(A 1 ...A 4 )- 2 Pof^ ) ..04A I 4)- 

As a consequence of this rescaling, and of the extra sum over A, there is a non-trivial amplitude associated to edges 
in the spin foam. Each edge contributes to the amplitude as 

where N\ to N4 are the colors of the four faces meeting at the given edge, and Ajvj,... ,at 4 is the dimension of the space 
of the interwiners between the representations Ni, . . . , A4. In conclusion, the amplitude of a Feynman diagram J is 
given by 

N f e v 

where the sum is over simple representations N of 50(4), and B denotes the Barrett-Crane vertex amplitude. 

Equivalently, as every edge connects two vertices, we can absorb the edge amplitude in the vertex amplitude and 
write Z in the standard form (]l|), where the vertex amplitude is 

n a 1/2 

A III "JV.1,... ,N i4 , „ / r ri\ 

v = 71 A ^2 OJV1...JV10) ( 59 J 

(AjVi ■ •• Ajvio) 

where Ni . . . Nio are the ten colors of the ten faces adjacent to the vertex v, and Nn . . . Nn, i = 1 ... 5 are the four 
colors of the four faces adjacent to i'th edge adjacent to the vertex v. 

We close this section with a comment. Unlikely the Barrett-Crane model, in the mode expansion of the model 
presented here, the field depends on five representations (four external and one intertwiner), which can be seen 
precisely as the quantum numbers of a "first quantized" geometry of a tetrahedron. 
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A. Bubbles in the spin foam formulation 



Now we are ready to show the finiteness of the 5-bubble amplitude (equation ( p5| ) in the previous section) directly 
in momentum space. As we mentioned before, the integration over p ao in ([35J) projects the amplitude into the trivial 
representations. The value of I5 is then given by the sum over colors of the amplitudes corresponding to the diagram 
on the right of Fig. [| in which the colors of the external faces are fixed to zero. Four of the colors vanish in each of 
the five vertices, and a typical vertex amplitude reduces to 



SooOOiVi...JV 6 



Ni N 2 N 3 
A 4 N 5 N 6 



(60) 



6j 



where the RHS denotes the 6 J— symbol defined in the appendix. Each of the internal faces is colored by irreducible 
representations iVi to iVio. The internal faces are triangles; therefore for each Ajv, contribution of the face, there will 
be 3 edge contributions ((A^ 2 ) 3 ). The factor in the numerator of ( |57| ) reduces to the general form Ao,Ni,n-,,n s = 1) 
since the space of intertwiners between 3 representations is one dimensional. Putting this together, we get from (pq) 



25= (A J v 1 ...A iVlD )(A JVl ...A J v 1D ) 

N!...N 10 



TT I N? N% N% 
111 Nl NV NS 



(61) 



Using the bound for the 6j-symbols given in ( A12| ), and noticing that each color appears in three of the vertices, we 
conclude: 

fa\< £ (A Afl ...A 7Vl0 )(A Ari ...A, Vl0 )- 6 (A Wl ...A Ari0 ) 3/2 = ]T (A Nl ...Ajv 10 )~' < 00. (62) 

iVi...iVia W1...JV10 

Next, consider the case in which the external colors are not zero. The finiteness of these terms follows from lemma 4 
in the appendix. The proof follows the similar s teps as before, simply replacing the bound for the 6j-symbols with the 
bound on the Barrett-Crane vertex amplitude (A24). More precisely, there are two kind of colorings corresponding to 
external and internal faces respectively. We denote by N e the colors labeling the ten external surfaces in the diagram 
of Fig. ||, while N l denotes the colors of the ten internal faces. Only internal colors are summed over. External 
colors label the surfaces shown on the left of Fig. 5. They appear in the propagator corresponding to the appropriate 
internal edge, thus contr ibuting with a factor A^ 2 to the amplitude. They also appear as arguments of two vertex 
amplitudes. Using (A24) this contribution will be less or equal than A NC . There will be no face contribution for 
external colors (the face contribution appears when one has a complete chain of propagators closing around a face 
so that all the S a p in ( |52] ) combine into S aa = An)- Therefore, external colors N e can appear in the bound for the 
amplitude only through the dimensionality of the space of intertwiners AN 1 ,... t N< L in each propagator associated to 
the ten internal edges. We denote the An 1 ,....n 4 , m as A ab for a,b = 1, . . . , 5 and a ^ b, according to Fig. ||. 

On the other hand, internal colors appear as face contributions (A N i)- There is also a contribution from the three 
edges of each face (Ajfi), and the contribution of the three corresponding vertices which are less or equal that A 3 ^ , 



according to ( |A24| ). 

From the considerations of the two previous paragraphs we obtain 



4 bc) (Ar...ivro) 



< 



Ni...Ni 



A 



N! 



, A 



If 



a ab 



(63) 



In order to construct a manifestly finite upper bound for the amplitude we need to obtain a bound for AN lt ...,N 4 - 
To find an upper bound for AjVi,...,JV4 we proceed as follows. The A's appear in (55) through 4> < N l N2N 3 N i Ai' wnere A 
labels the elements of an orthonormal base of intertwiners between the four representations N\ . . . N±. Assume that 
Ni is less or equal than the other N. We chose a basis in which for example Ni and N2 are in the same 3-intertwiner 
QNtNzA in (Ta|) g ince jy i5 jy 2) and A have to be s0 ^ compatible it follows that N 2 - Ni < A < Ni + N 2 . There 

satisfies Ajv^ 



are 2Ni + 1 = yAjVi A's who satisfy this condition. However, there are additional compatibility conditions on A 
so we conclude that the number of possible A's denoted by A^ 
smallest N we can also write 



,N 4 



< ^/Ajvi, since Ni was the 



A 



iVi,...,JV 4 



(64) 



for i = l,2, 3, 4. The dimensionality of the space of intertwiners is independent of the chosen basis and the previous 
inequality has a basis independent meaning. We can use the colors of the ten internal faces N l to bound the ten A ab 
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in the previous equation. Finally, according to ( p4j ) we obtain 

_ 3^ 

< ( A ^ • • ■ A Ni ) ' <M - ( 65 ) 

The same kind of inequality can be derived for A\ c . 

V. CONCLUSION 

We have presented in this paper a spin foam model, possibly related to Euclidean quantum gravity, in which the 
fundamental amplitudes, divergent in the BF and Barrett-Crane models, are finite. The model is obtained as a 
modification of the interaction term of quantum BF theory, formulated on the group manifold. The modification is 
an implementation, a la Barrett-Crane, of the constraint that reduced BF to general relativity. 

We have not proven finiteness of all amplitudes at all orders in the expansion in 2-complexes, but we suspect 
that by using the methods introduced here the proof could be possible. For each particular amplitude, it is easy 
to construct a finite bound of the amplitude, using the inequalities given in the paper. The difficulty consists of 
finding a general way of bounding the degeneracy of the vertex amplitude given by the dimensionality of the space of 
intertwiners associated to the 5 edges converging at each vertex. In addition, the inequalities given in the appendix 
might be used to study the convergence of the full Feynman series, or sum over 2-complex. 

It would be crucial to get a better understanding of the classical limit of the model. In particular, it would be 
interesting to know whether the amplitude of a spin foam approaches the exponential of the Einstein-Hilbert action 
of the four-geometry that the spin foam approximate. 

More in general, we think that the techniques and methods introduced here could be useful for analyzing divergences 
in general spin foam models. 



aW(nz...n; ) 
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APPENDIX A: IRREDUCIBLE REPRESENTATIONS OF SO (4) 

In this appendix we review some properties and definitions in the theory of irreducible representations of SO (4) 



(on this see 31 ), and we state and prove some i neq ualities used in the paper to show the finiteness of amplitudes in 



previous sections. We follow the conventions of |E3[. 

Given g e SO (4) we denote by D^{g) the representation matrix corresponding to the irreducible representation 
of order A. Integration over SO(4) or the SO(3) subgroup H is performed with the normalized Haar measure of the 
group and the subgroup respectively. The integration of two representation matrices is given by 



I 



dg D<$(g)D%\(g) = -Ls AA ' <W <W, (Al) 

SO(4) 

where A(N) denotes the dimension of the representation. In the case of SO(4) we can choose a basis in which 
matrices are orthogonal, and the bar can be dropped from the previous equation. The integral of the product of 
three group elements is 

/ * D™ (g)Dl N J 2 (g)D^ (g) = C™. (A2) 

SO(4) 

Here C£j^~ A are normalized intertwiners (Wigner 3-j symbols) between three representations of SO (4); that is 
Caiaza ^aia 2 a = 1 • The intertwiner from the tensor product of two representations Ni , N% to a representation 
A, if it exists is unique. 
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The intertwiners between four representations are of great importance in our calculation. However, in general they 
are not unique and rather form a vector space. An orthonormal base can be defined as follows: 



A = C%%\ (A3) 

With these definitions equation (|A^ ) generalizes to the case of the integration of four representation matrices to 



/ ^<^(.9)...<£( 5 ) = E C ^: K 



N 4 A 



(A4) 



SO(4) 



Another important equation corresponds to the integration of one representation matrix over a sub-group SO(3) C 
50(4), namely 



dhDWh)=wWv>™ 



(A5) 



H=SO(3) 



where represents the unit vector in the irreducible representation of order N left invariant by the action of the 
subgroup H (iVa is non vanishing only if TV is simple). Equation (A5) defines the projector into that one-dimensional 
vector space. 

As mentioned, invariant vectors exist only in simple representations. As a consequence the projection of the 
intertwiner C^ 1 '^ 4 w w 7l . . . w l4 vanishes unless all the iVj and N are simple. In this case its value is given by 



(-1NL..N4N 
°7i-74 w ii ■■■ w n 



1 



y /A(N 1 )...A(N A y 
Finally we give the definition of the Barrett-Crane intertwiner: 



71—74 ~ / j 7i ... 1 



N 4 N 

74 



(A6) 



(A7) 



The previous is the un-normalized Barrett-Crane intertwiner as originally defined in [Q which is shown to be unique 
up to scaling in |p3|,p2[ . 



1. Some lemmas 



Now we state and prove the lemma that plays an important role in showing the finiteness of the 1-bubble amplitude 
in the modified Barrett-Crane model. We define the following distribution 

Z>(ffi,0>,Sa)= / ^dft^/3i^i/?r 1 /?2^2/3 2 ~ 1 5i/33^/3 3 "V4/i4^ 1 52/55^5/?5" 1 /36^/3 6 " 1 53) (A8) 

Lemma 1: 

The distribution T>(gi, 32733) is a bounded function over 5'0(4) 3 . 



(A9) 



Proof: In order to prove the lemma we expand (|A8| ) using 

%) = £A(AOTr(£W( 5 )). 



N 



We integrate over /3j and hi and make use of equations ( |A5[ ) and (Al). The mode expansion of 2?(<?i, g 2 , <? 3 ) results: 

^(51,52,33)= E A^A) Tr(jD(A) (3l52ff3)) - (A10) 



Simple A 



Now 
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Simple A ^ ' 



< E a^a)|t^ (a) (^3))|< £ A^y< 

Simple A Simple A 



OO, 



(All) 
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where we have used that in our orthogonal representation \Tr(D A (g))\ < Aa. Therefore, T>(gi, g 2l g 3 ) is bounded on 
50 3 (4) □ 

By means of the previous lemma we were able to find a bound to the 1-bubble amplitude directly in configuration 
space. To show the finiteness of the 5-bubble amplitude in the previous section, we made use of the mode expansion 
of the amplitude and the following lemma|. 
Lemma 2: 

The 50(4) 6j-symbols defined in terms of normalized intertwiners (normalized Wigner 3j-symbols C^a^al 3 °^ 
(Ia5)) satisfy the following inequality: 
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(A12) 



where A at denotes the dimension of the irreducible representation of order N. 
Proof: The 6j-symbol is defined in terms of C^a^ 3 as 



Nl N 2 N3 1 _ r ,N 1 N 2 N ar ,N 1 N lL N 6r iN2N lL N &r ,N 3 N 6 N 5 

N 4 N 5 N 6 J ~ ' 
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where summation over repeated indices is understood. The assertion of the lemma is proven by means of calculating 
the following integral in two different ways. 



S 



dgidtodgadgifa^DfafifaDfa)^) (D(g 2 )^D(g 2 )^D(g 2 )2 

(i?(93)^^(53)^ 4 i5( 93 )^ 5 )(i?( 9 4)^ 3 i)( 5 4)^ 6 i?( 9 4)^ 

We can rewrite the previous equation using the representation property D(g)D(f) = D(gf) as 

S = J d 9i Tr [D^(g l9 ^)} Tr [D»* fag^ 1 )] Tr [D»» fagf 1 )] 

Tr [D^(g 2 g^)] Tr [D»< (g 3 g^)] Tr [D^(g 2 g^)] . 



(A14) 



(A15) 



The fact that we are using orthogonal irreducible representations of 50(4) implies that |Tr[D Ar (<7)] | < A at- Combin- 
ing this bound for the trace with the normalization of the SO (A) Haar measure in the previous equation we conclude 
that 



|5| < A Nl A N2 A Ns A Ni A Ns A Nfj . 



On the other hand, using equation (A2), and the definition of the 6j'-symbol above we obtain 

2 



S = 



N x N 2 N 3 
Ni N 5 N 6 



(A16) 



(A17) 



6j 



which concludes the proof □ 
Lemma 3: 

The 50(4) 15j-symbols defined in terms of normalized intertwiners (C^ ( ^ 2 2 ( ^^ i ' A of (A4)) satisfy the following 
inequality: 
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Nq n 7 n$ n 9 n w 
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(A18) 



15j 



where 



We have not found a reference to this lemma in the literature, although we suspect is must be a known result due to its 
simplicity. 
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(A19) 



This lemma generalizes the previous one; however for simplicity we have proven explicitly the first. The proof 
of the current lemma follows the analogous path as the previous one with some additional indices. One starts by 
defining an integral analogous to S in lemma(2); namely, 



S = J ^i^3^5^Si)^ 

D(92)^D(g 2 )^D(g 3 )^D(g^ 
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(A20) 



where now one has an integration over five group variables of the appropriate product of twenty represent atio n 
matrices. One can relate the value of the integral to the sum of 15j-symbols squared over the A's by means of (A4). 
Finally, one finds a bound to the integral using the fundamental inequality | Tr[i? JV '(g)] | < A^r- 
Now we state two corollaries of the previous lemma. First, we have that 



Ni N 2 N 3 N 4 N 5 
N 6 N 7 N s N 9 N 10 
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and second, from the definition of the Barrett-Crane vertex we conclude that 



\Bn u 
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An important consequence of (|A20| ) is that it can also be used to directly find a bound for the Barr ett-C rane vertex. 
Take the integral defined in flA20|) in which the $ are not contracted between them self as in (|A20|) but rather 
contracted to the appropriate set of normalized invariant vectors w^* in each representation. Using (|A4|) and ( |A6| ) 
we obtain, on the one hand 
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Nt,...,N 10 
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(A23) 



However, on the other hand, the absolute value of integrand contracted with all the w's is less or equal than one 
since it can be written as a product of terms of the form \w li D{g)^w^\ < 1. Therefore, we have proven the following- 
lemma: 

Lemma 4: 

The Barrett-Crane vertex amplitude satisfies the following inequality: 
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